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ABSTRACT 


Many  current  army  combat  simulation  models  require 
digitized  terrain  inputs  for  their  execution,  resulting 
in  a costly  and  time-consuming  operation.  A methodology 
to  generate  parametric  terrain  as  a function  of  the 
number  of  hill  masses  and  the  height,  location,  spread, 
and  slope  of  the  hills  is  presented  in  this  paper.  Modi- 
fied bivariate  normal  (MBVN)  distribution  function  and  a 
linear  function  cutting  the  MBVN  function  are  employed  to 
achieve  the  terrain  representation.  Analysis  of  various 
terrain  parameter  values  was  conducted  and  the  significant 
trends  resulting  from  the  study  are  described.  Finally, 
recommendations  for  future  enhancement  of  the  model  are 
given. 
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I . INTRODUCTION 


Many  current  army  combat  simulation  models  require 
digitized  terrain  inputs  for  their  execution.  This  process 
is  very  costly  and  requires  a significant  amount  of  computer 
storage.  In  addition,  computation  of  line-of-sight  and 
mobility  require  a great  deal  of  computer  running  time. 

The  current  methodology  utilized  to  represent  terrain 
is  also  questionable  from  a statistical  point  of  view. 

A Vector  Research  Incorporated  (VRI)  study  on  terrain  line 
of  sight,  concluded  that  "...  present  and  past  army  study 
results,  based  on  the  analysis  of  combat  results  on  a very 
limited  sample  of  terrains,  may  have  been  determined  by 
the  terrain  selection  process  and  not  by  the  actual  weapon 
system  or  force  design  differences."  More  specifically 
it  stated  the  following; 

A.  "There  is  extreme  sensitivity  in  combat  model 
results  as  the  scenarios  (terrain  and  movement  assumption) 
are  varied,  even  when  variation  is  within  a class  of 
scenarios  chosen  for  their  a priori  equivalence." 

B.  "This  sensitivity  can  be  slightly  reduced,  but 
remains  extreme  (with  probabilities  of  win  estimable  only 
within  plus  or  minus  25%)  even  when  battle  results  are 
used  to  redesign  scenarios."^ 


■‘'Farrell,  Robert  L.  and  Richard  J.  Freedman,  Investigation 
of  the  Variation  of  Combat  Model  Predictions  with  Terrain  Line 


_t.  Vector  Research  Incorporated,  19  7 
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These  results  imply  that  sufficient  replications  of  each 
type  of  terrain  should  be  run  in  order  to  reach  a satisfactory 
statistical  level  of  significance.  For  VRI's  analysis  at 
least  fifty  replications  of  each  type  of  terrain  were  used. 

In  view  of  the  problems  stated  cibove , research  was  xinder- 

2 

taken  by  Major  Christopher  Needels  to  develop  a new  method- 
ology for  the  representation  of  terrain.  In  particular,  a 
modified  bivariate  normal  distribution  function  was  utilized 
to  generate  hill  masses  in  a wide  variety  of  configurations. 
Line  of  sight  and  movement  routines  were  also  developed  to 
evaluate  various  intervisibility  parameters  over  many  terrain 
configurations  for  very  low  cost. 

This  thesis  represents  an  extension  of  the  models 
developed  by  Needels  to  generalize  the  terrain  generation 
routines.  In  particular,  methodology  was  developed  to 
generalize  the  hill  mass  generation  from  symmetric  hill 
shapes  to  non-symmetric  configurations.  The  methodology 
developed  for  this  generalization  is  described  in  Chapter  III. 

In  order  to  evaluate  the  effects  of  the  number  of  hills, 
hill  height,  and  slope  on  various  intervisibility  parameters, 
an  experiment  was  designed  to  measure  these  effects.  A 
description  of  the  design  methodology  and  results  of  the 
analysis  are  given  in  Chapter  IV.  Finally,  the  conclusions 
drawn  from  the  current  analysis  cind  recommendations  for 
future  research  are  presented  in  Chapter  V. 

2 

Needels,  Christopher  J. , Parameterization  of  Terrain  in 
Army  Combat  Analysis,  Mas  ter ' s Thesis,  Naval  Postgraduate 
School,  Monterey,  California,  1976. 
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II.  REPRESENTATION  OF  TERRAIN 


A.  HIGH  RESOLUTION  COMBAT  MODEL  REPRESENTATION 

The  classic  approach  to  representing  terrain  in  land 
combat  models  is  to  utilize  digitized  terrain  evaluation 
data  from  Waterway  Experiment  Station  (WES)  for  the 
particular  area  of  interest. 

The  Dynamic  Tactical  Simulator  (DYNTACS)  model  serves 
as  a good  example  of  the  state  of  the  art  in  terrain 
representation  and  use.  DYNTACS  is  a two-sided  dynamic 
model  of  battalion  level  combat  representing  the  detailed 
interaction  of  elements  in  a combined  arms  environment. 

DYNTACS  utilizes  100  meter  grid  squares  for  which 
macro-terrain  elevations  are  determined  from  WES  tapes. 

The  simulation  divides  each  grid  square  diagonally,  thus 
providing  a series  of  adjoining  triangular  terrain.  The 
entire  battlefield,  therefore,  is  represented  as  a surface 
of  equal  size  triangles  which  vary  in  slope  depending  on 
the  elevation  at  their  corners. 

Needels  (pp.  9-10)  provides  a concise  description  of 
the  DYNTACS  methodology  for  line  of  sight  calculation  as 
follows . 


From  the  macro-terrain  data,  line  of  sight 
between  any  two  opposing  elements  is  computed. 
This  is  accomplished  by  first  computing  the 
angle  between  the  horizontal  and  a straight 
line  drawn  between  an  observer  and  target 
(0-T  Line) . The  program  then  conducts  a 
search  of  the  terrain  along  the  path  of  the 
0-T  line  to  see  if  any  macro- terrain  is 
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higher  than  the  0-T  line  itself.  This  is 
accomplished  by  comparing  the  angle  of  the 
0-T  line  with  the  angle  eibove  horizontal 
of  the  Observer-Terrain  line.  If  the  latter 
angle  is  larger,  there  is  no  intervisibility. 
Over  the  duration  of  a battle  with  numerous 
elements  this  calculation  may  be  made 
thousands  of  times.  Consequently,  not  only 
is  the  time  to  prepare  the  tapes  high,  but 
also  the  time  to  compute  lines  of  sight  once 
the  terrain  is  input  to  the  model. 
Considereible  effort  by  developers  and  users 
of  this  model  has  been  expended  to  stream- 
line this  subroutine. 


B.  SYMMETRIC  REPRESENTATION 

3 

The  basic  models  developed  by  Needels  (hereafter  referred 
to  as  SIMTER)  which  provide  the  basic  foundation  for  the 
methodology  developed  in  this  thesis  are  described  in  this 
section.  The  basic  motivation  for  SIMTER  was  to  develop  a 
random  terrain  representation  which  could  be  replicated 
quickly  but  at  the  same  time  be  representative  of  a particular 
"type”  of  terrain.  The  approach  adopted  for  SIMTER  generated 
a variable  number  of  hill  masses.  The  shape  and  height  of 
the  hills  could  be  varied,  but  each  hill  was  symmetric,  since 
it  was  generated  directly  from  the  modified  bivariate  normal 
(MBVN)  density  function. 

SIMTER  has  two  parts.  The  first  part  is  the  main 
progreun  which  creates  terrain  and,  at  the  discretion  of  the 
user,  plots  both  a three-dimensional  drawing  and  a contour 
map.  Used  strictly  for  terrain  generation,  it  can  be  used 
as  a preprocessor  for  other  combat  models  by  producing  grid 
points  and  elevations  similar  to  those  provided  on  computer 

^Needels , Parameterization  of  Terrain  in  Army  Combat 
Analysis . 
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tapes  by  WES.  The  second  part  of  SIMTER  is  a movement/line 
of  sight  routine  which  moves  a target  along  specified  routes 
across  the  generated  terrain. 

C.  A NEW  REPRESENTATION  METHODOLOGY 

The  methodology  developed  in  this  thesis  provides  for 
the  generation  of  hi21  masses  utilizing  MBVN  density 
function  to  represent  masses  that  are  not  symmetric.  This 
development  provides  a significant  increase  in  flexibility 
to  represent  hill  masses  more  realistically.  The  mathematical 
development  of  this  methodology  is  presented  in  the  next 
chapter,  after  which  actual  results  from  executing  the  model 
are  given. 
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III.  THE  EXTENDED  MODEL 


A.  UNSYMMETRIC  TERRAIN 

The  unsymmetric  slope  of  terrain  can  be  generated  by 
a cutting  method  using  the  MBVN  as  follows: 


Z 


where 


i 

I 


MH • EXP  \ - 


X(l-p^) 


[(- 


x-y 


X.  2 


2 ( 


MH:  maximum  height  of  hill 

Z : height  of  hill 

y^,y^:  spread  of  hill 
p : ellipse  factor  of  hill 


The  unsymmetric  slope  of  the  hill  can  be  represented  by  the 
vertical  value  of  the  MBVN  subtracted  from  the  vertical  value 
of  a linear  function  as  shown  in  Figure  1 by  the  value  HI. 


FIGURE  1.  UNSYMMETRIC  SLOPE  GENERATION 
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Two  basic  questions  must  be  addressed: 

1.  Where  do  we  start  the  cutting  of  the  MBVN  function? 

2.  How  do  we  compute  the  initial  cutting  position? 

If  M is  the  direction  of  the  cutting  line  (Line  A) , then 

the  cutting  plane  is  through  the  location  of  the  center  of 

the  hill  (y  ,y  ) as  shown  in  Figure  2. 

X y 


FIGURE  2.  CUTTING  LINE  DETERMINATION 

If  the  initial  point  of  cutting  is  (XIN,YIN)  and  HI  is 
height  at  the  initial  point,  this  point  is  on  the  cutting 
plane  and  also  on  the  MBVN  distribution  function.  Thus,  if 
we  consider  the  graph  in  two  dimensions  as  shown  in  Figure  2, 
Line  B is  tangent  to  the  circle  at  the  point  that  the  vertical 
value  on  the  MBVN  function  is  HI  and  the  slope  of  Line  B 
is  -1/M.  If  the  initial  point  of  the  cutting  is  (XIN,YIN) 
this  point  is  on  Line  B,  Line  A cind  the  circle  at  a height 
HI  on  the  MBVN  function.  We  can  compute  the  initial  point 
(XIN,YIN)  as  follows. 


13 


1.  0=0  Case 


1 YIN-u 

- i[(— ~2i)2  + (__J:)2. 


HI  = MH*EXP<-  i[( 

/ 2 o,  Oy 


XIN-y  _ YIN-y  _ 

-2  ln(HI/MH)  = ( + ( 1)2 

°x 


From  (2),  the  derivative  is  computed  as  follows; 


dYIN 

dXIN 


XIN-y 
( 7^) 


YIN-y 
( 7-) 


From  (2)  and  (3) 
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* * / 1 / _ /-  \ ^ 
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y y I 


-2  ln(HI/MH) 
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If  the  hill  is  cut  from  left  to  right  and  ° t then 


Thus,  if  the  direction  of  cutting  and  HI  are  selected,  then 
the  starting  point  can  be  computed.  Care  must  be  taken  in 
the  selection  of  the  value  of  HI.  If  HI  is  too  small,  then 
the  hill  will  have  an  unrealistic  shape.  The  value  of  HI 
affects  the  shape  of  the  hill  as  shown  in  Figure  3. 


FIGURE  3.  EFFECT  OF  HI  ON  THE  HILL  SHAPE 


15 


2 . 0 < p < 1 Case 
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From  (13)  the  derivative  is  computed  as  follows: 
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From  (14) 
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Let  C 


From  (13)  , 


-2(1  - 
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Then 


XIN  = 
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X X 
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(15)  and  (16) , 


p ) ln(HI/MH)  = (1 
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Next,  the  height  of  intersecting  hill  masses  must 


be  determined  as  illustrated  in  Figure  5. 


XA 


FIGURE  5.  INTERSECTING  HILL  MASSES 

HB  defines  height  of  hill  "B"  after  cutting. 

HA  defines  height  of  hill  "A"  after  cutting. 

Height  of  hill  "A"  is  greater  than  that  of  hill  "B" 
at  the  point  XA  but  height  of  hill  "B"  is  greater  than  that 
of  hill  "A"  after  cutting.  If  there  is  more  than  one  hill 
located  close  to  each  other,  those  hills  must  be  adjusted 
by  comparing  the  height  after  cutting  each  one . The 
computation  of  the  height  of  the  terrain  at  any  point  is 
described  in  the  following  section. 

B.  TERRAIN  HEIGHT 

The  height  of  the  hill  which  is  generated  by  the  MBVN 
distribution  corresponds  to  the  height  of  a hill  which  is 
generated  by  the  cutting  method.  The  height  of  a hill 
at  euiy  point  is  computed  as  follows: 
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where  0 = cutting  angle  of  the  x-axis 

(})  = cutting  amgle  of  the  y-axis. 

The  height  of  the  hill  after  cutting  is 

if  Z < HI  and  ZCUT  < 0 

if  Z < HI  and  ZCUT  > 0 

if  Z > HI  and  ZCUT  > 0 

if  Z > HI  and  ZCUT  > 0 

and  Z > ZCUT+HI 

if  Z > HI  and  ZCUT  > 0 

and  Z < ZCUT+HI 

If  more  than  one  hill  exists  at  a point  on  the  terrain,  the 
greatest  height  after  cutting  is  used. 

C.  MOVEMENT  LINE  OF  SIGHT 

The  procedure  of  computing  the  line-of-sight  between 
any  two  points  on  the  terrain  is  described  in  this  section. 
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If  the  position  of  the  observer  is  (OX,OY)  and  the  position 
of  the  target  is  (X,Y) , then  the  slope  of  the  Observer- 
Target  (0-T)  Line  is  given  by 


ML  = 


X - OX 
Y - OY 


(28) 


If  the  intersection  point  of  a hill  and  the  0-T  Line  is 
(XI, YI),  then  HI  eind  the  0-T  Line  are  perpendicular  as 
shown  in  Figure  6 . 


0^ 

(DX, oy) 


cx  ,> ) 


Thus , 


FIGURE  6.  INTERSECTION  OF  0-T  LINE 
WITH  HILL  MASS 


XI  = 


YI  = 


ML  •OX-ML'OY  + y +ML-y 

_2S Y. 

1 + ML^ 

ML  (XI  - OX)  + OY 


(29) 


(30) 
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If  there  is  more  than  one  hill  on  the  0-T  Line,  the 
dominant  hill  (after  cutting)  is  used  to  compute  the  percent 
of  the  target  covered  relative  to  the  observer.  Target 
movement  along  the  surface  slope  of  the  hill  is  computed 
as  follows: 


ZH  = MH*EXP  {- 


2(l-p^)  -X 


-)^  - 2p( -)  { il) 

a ^ 0 a 

XV 


y"^v  2 
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(33) 
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If  DEL  is  the  direction  of  target  movement,  then 


DEL 


3ZH 

9x 


cos  r 


+ 


9ZH 

3y 


sin 


r 


(34) 


If  the  speed  of  the  target  is  constant  "V"  on  the  surface, 
then 


velocity  = ■ (35) 

y 1 + (DEL)^ 

If  the  speed  of  the  target  is  dependent  on  the  slope  of 
the  surface,  then 

velocity  = V + W(DEL)  (36) 

where  W is  the  speed  factor  of  slope. 
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IV.  TEST  AND  EVALUATION  OF  MODEL  RESULTS 

A.  DESIGN  AND  CUTTING  ANGLE  METHODOLOGY 

The  methodology  for  parametrically  describing  terrain 
as  a function  of  the  number  of  hills,  the  center  and  spread 

of  each  hill,  and  the  slope  was  developed  in  the  previous 

chapters.  In  this  chapter  the  effect  of  the  number  of  hills, 
height  of  hills,  and  the  slope  on  various  intervisibility 
measures  is  investigated.  Three  random  terrains  were 
generated  (with  each  terrain  having  six,  twelve,  and  eighteen 
hills)  resulting  in  nine  basic  terrain  configurations. 

For  each  configuration  three  average  hill  heights  (100, 

200,  and  300  meters)  were  utilized. 

Finally,  for  each  of  the  27  resulting  configurations, 
four  slope  cutting  methods  were  used  as  follows: 

1.  Down  slope-cut  from  observer  side  to  target  side. 

2.  Side  slope-cut  from  right  to  left  as  viewed  from 

the  target. 

3.  Up  slope  — cut  from  target  side  to  observer  side. 

4.  Both  slope-symmetric  slope  as  viewed  from  any 
position . 

A total  of  108  vinique  terrain  configurations  were 
generated  for  the  investigation.  For  each  of  the  config- 
urations, the  mean  and  standard  deviation  describing  the 
center  and  spread  of  each  hill  is  specified.  In  addition, 
the  cutting  angle  for  determining  slope  must  also  be 
specified. 
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If  the  cutting  angle  is  constant,  then  the  hill  can  be 
described  as  follows: 

HILL  = F(Peak,  Position,  Standard  Deviation,  Slope) 

In  the  constant  cutting  angle  case,  only  the  slope  is 

variable,  with  peak,  standard  deviation,  and  position  < 

constant.  In  other  words,  if  the  slope  varies  in  Figure  7, 
the  standard  deviation  of  the  resulting  distribution  also 
varies . 


FIGURE  7.  CUTTING  ANGLE  AND  SLOPE 


The  objective  of  the  methodology  is  to  maintain  a constant 
standard  deviation  of  the  distribuition  function  resulting 
from  each  cutting  angle,  so  that  only  the  slope  of  the  hill 
varies.  Otherwise  the  determination  of  how  slope  affects 
the  intervisibility  parameters  would  not  be  possible,  since 
it  would  be  confound  with  the  standard  deviation  (spread) 
of  the  hills. 
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The  problem,  therefore,  was  to  determine  a method  to 
use  various  cutting  angles  and  maintain  a constant  standard 
deviation.  After  extensive  investigations  into  various 
approaches  to  this  problem,  no  method  was  found  which  would 
mathematically  guarantee  a constant  standard  deviation. 

The  method  shown  in  Figure  8 was  adapted  as  the  best 
approximation  method  available. 


figure  8.  CUTTING  ANGLE  METHODOLOGY 

In  Figure  8,  HT  is  the  maximum  height  of  the  hill  after 
cutting  with  Line  A.  Note  that  the  cutting  angle  is  now 
specified  in  reference  to  the  center  line  of  the  original 
bivariate  normal  hill  mass  distribution.  In  order  to 
generate  a symmetric  hill  mass  with  approximately  the  same 
standard  deviation  (Line  C)  the  cut  is  made  maintaining 
a constant  maximum  height,  HT,  as  in  Figure  8. 

Several  test  runs  made  using  this  methodology  indicated 
that  the  resulting  standard  deviations  were  very  nearly 
equal  for  cutting  erngles.  However,  veiry  little  change  in 
slope  results  for  cutting  angles  greater  than  twenty  degrees. 
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B.  OBSERVER- TARGET  POSITIONS 


For  each  terrain  an  observer  position  and  initial  target 
position  were  selected.  Two  routes  from  the  target  to  the 
observer  position  were  selected  in  accordance  with  accepted 
tactics  of  advance  routes.  Selected  terrains  indicating 
the  observer-target  position  and  selected  routes  are  shown 
in  Figures  12  through  17.  in  order  to  assure  that  the 
"number  of  hills"  factor  was  properly  cinalyzed,  the  same 
routes  of  approach  on  each  generated  terrain  were  used, 
independent  of  the  number  of  hills. 

C . TARGET  MOVEMENT 

Two  cases  of  target  movement  were  considered  in  the 
analysis . 

1 . Constant  Speed 

The  target  moves  along  the  route  at  a constant  speed 
of  six  meters  per  second,  independent  of  hill  slopes.  The 
constant  speed  runs  serve  as  the  base  case  for  measurement 
of  the  percent  of  time  intervisible  along  the  route. 

2 . VaricJale  Speed 

The  speed  of  the  target  is  a function  of  the  terrain 
slope  being  traversed.  The  speed  of  the  vehicle  is  computed 
by  the  relationship  as  follows; 

velocity  = 6.0  + 4.0 (slope) 

This  slope  is  the  directional  derivative  value  at  the 
instantaneous  target  position  on  the  hill. 
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D.  TARGET  INTERVISIBILITY 


Five  dependent  variables  describing  target  intervisibility 
are  computed  for  each  advance  route  as  follows : 

1.  Percent  of  time  intervisible  along  the  route, 

2.  The  number  of  intervisibility  segments  along  the  route, 

3.  The  average  intervisibility  segment  length, 

4.  The  average  intervisibility  segment  time,  and 

5.  The  total  length  of  intervisibility  over  the  route. 
Intervisibility  pareuneters  are  computed  by  utilizing  the 
line  of  sight  algorithm  at  fixed  time  intervals,  TL,  during 
the  movement  run.  For  this  analysis,  TL  was  set  equal  to 
one  second.  A sample  output  for  one  compute  movement  trace 
is  given  in  Figure  9. 

E.  EXPERIMENTAL  DESIGN 

The  objective  of  the  analysis  is  to  determine  the  effects 
of  the  number  of  hills  (A) , the  height  of  hills  (B)  and  the 
slope  of  the  hills  (R)  on  target  intervisibility.  To 
accomplish  this  objective,  an  experiment  was  designed  as 
described  by  the  following  design  matrix. 

SLOPE 

# of  Hills  Height  Down  Side  Up  Both 

100 

6 200 

300 

100 

12  200 

300 

100 

18  200 

300 
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FIGURE  9.  OUTPUT  OF  DATA 
Definition  of  Varia±>les 

N Number  of  hills 

XMEL,YMEL  Spread  of  hills 

HITE  Maximum  height  of  hill 

RUFX, RUFY, RUFP  Different  amount  of  the  spread  and  height 
among  hills 

INCX,INCY  Grid  increment  for  the  terrain  matrix 

XMAP,YMAP  Size  of  the  grid  system  in  meters 

XGRID,YGRID  Center  position  of  hill 

RIGHT  Maximum  height  of  hill  in  modified  bivariate 

normal  distribution 

XSTD/YSTD  Spread  of  hill  (standard  deviation) 

XSLOP/YSLOP  Cutting  angle 

DIRECTION  Cutting  direction 

XINIT/YINIT  Initial  position  of  cutting  plane 

HEIGHT  Hieght  of  hill  at  the  center  position  on  MBVN 

after  cutting 

RX,RY  Coordinates  of  approaching  route  of  target 

TIME  Movement  time 

MAP  COORDINATES  Target  position  at  given  time 

ELEVATION  Target  elevation  on  the  terrain 

SPEED  Velocity  on  the  terrain  at  the  given  time 

VISIBLE-  1 means  visible  cmd  0 means  invisible 

HLIN  The  different  height  of  0-T  line  at  the 

maximum  height  of  surface  of  the  terrain  on 
the  0-T  line,  and  the  lower  height  between 
observer  and  target 
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DZ 


Maximum  height  of  surface  of  the  terrain  on  the 
0-T  line 


KL  The  hill  which  is  maximiam  height  between  the 

observer  emd  target 

DLIN  Difference  maximum  height  of  surface  between 

observer  emd  target,  and  lower  height  between 
observer  and  target. 
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FIGURE  9 CONTINUED 
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The  appropriate  linear  model  is  as  follows: 


where 


„ f t B.  t + ABi.  + AR.^  + BR.^  * ABR. 


Y 

ijk£ 

y 

A 


intervisibility  dependent  variable  measure 
grand  mean  of  the  population 

nvimber  of  hills,  where 

1-6  hills 
i = 2-12  hills 

3-18  hills 


B = average  height  of  hills,  where 

1-100  meters 
j = 2-200  meters 

3- 300  meters 

R = direction  of  slope  of  hills,  where 

1-down  slope 
, _ 2-side  slope 
^ ~ 3-up  slope 

4- both  (symmetric) 

Within  each  cell  of  the  matrix,  there  are  six  replications 


consisting  of  two  routes  on  each  of  the  three  basic  terrain 


configurations . 


For  the  first  treatment  levels  of  slope  defined  above, 
a cutting  angle  of  fifteen  degrees  was  used,  the  symmetric 
case  was  generated  with  a cutting  euigle  of  zero  degrees 
using  the  methodology  previously  described.  Four  analyses 
were  conducted  as  representative  of  the  large  number  of 
possible  sitatuations . The  measurements  selected  were 

1.  Percent  of  time  intervisible  for  constant  speed 
and  for  variable  speed,  and 

2.  Number  of  intervisibility  segments  for  constant  ^peed 
and  for  variable  speed. 
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The  results  for  each  selected  analysis  are  described  in 
the  following  section. 


F.  ANALYSIS  OF  RESULTS 

1 . Percentage  of  Intervisibility  Time 
a.  Constcuit  Speed 


1/8 

(x  ' Treinsgeneration  Data) 

SUMMARY  TABLE  - ANOVA 


Source 

DF 

ss 

MS 

F Ratio 

# of  hills 

2 

0.0437 

0.0218 

23.6417 

A 

height  of  hill 

2 

0.0117 

0.0058 

6.3209 

B 

slope  of  hill 

3 

0.0387 

0.0129 

13.962 

R 

AB  interaction 

4 

0.0006 

0.0081 

0.1527 

AR  interaction 

6 

0.0026 

0.0004 

0.4759 

BR  interaction 

6 

0.0016 

0.00027 

0.2943 

ABR  interaction 

12 

0.0012 

0.0001 

0.1114 

Error 

179 

0.1655 

0.0009 

Each  of  the  three  main  effects 

above  are 

significant 

: the  0.01  level 

, but 

no  interaction  terms 

are  significant. 

The  following 

teible 

gives  the 

mean  and  standard  deviation 

for  each  treatment  level  of  slope. 
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Slope 

Mean 

Std.  Dev 

down 

0.6130 

0.1539 

side 

0.5214 

0.1701 

up 

0.6926 

0.1629 

both 

0.5573 

0.1595 

Note  that  the  up  slope  has  the  greatest  percentage  of  inter- 
visibility time;  down  slope,  both,  amd  side  slope  following 
in  that  order. 

The  following  tadsle  gives  the  mean  and  staindard 
deviation  for  the  treatment  levels  of  number  of  hills  and 
height  of  hills  for  the  constant  speed  case. 


# of  Hills 

Mean 

Standard  Deviation 

6 

0.6976 

0.1620 

12 

0.5462 

0.1620 

18 

0.5272 

0.1490 

Height  of  Hill 

Mean 

Standard  Deviation 

100 

0.6475 

0.1775 

200 

0.5844 

0.1684 

300 

0.5556 

0.1628 
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b.  Speed  which  Depends  on  Slope 


Transgeneration  Data) 


SUMMARY  TABLE  - ANOVA 


Source 

DF 

SS 

MS 

F Ratio 

# of  hills 

2 

0.5152 

0.2576 

34.286 

A 

Height  of  hill 

2 

0.2041 

0.1020 

13.582 

B 

Slope  of  hill 

3 

0.3830 

0.1277 

16.991 

R 

AB  interaction 

4 

0.0033 

0.0008 

0.1090 

AR  interaction 

6 

0.0381 

0.0064 

0,8450 

BR  interaction 

6 

0.0165 

0.0028 

0.3670 

ABR  interaction 

12 

0.0091 

0.0008 

0.1000 

Error 

179 

1.3449 

0.0075 

Each  of  the  three  main  effects  cd)ove  are  significant  at  the 
0.01  level,  but  no  interaction  terms  are  significant. 

c.  Constant  and  Dependent  Speed. 

The  ANOVA  for  comparing  constant  and  dependent 
speed  is  given  in  Figure  10.  The  treatment  level  data  is 
as  follows; 

Meam  Standard  Deviation 

Constant  Speed  0.569  0.1732 

Dependent  Speed  0.5393  0.1921 
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^2C  CASES  WERE  PROCESSED. 

108  CASES  I 25.1  PC7J  WERE  MISSING 


2 . Number  of  Intervisibility  Segments 
a.  Constant  Speed 


SUMMARY  TABLE  - ANOVA 

1/2 

(X  Transgeneration  Data) 


Source 

DF 

SS 

MS 

F-Ratio 

# of  hill 

2 

8.3415 

4.1708 

82.281 

A 

Height  of  hill 

2 

0.2402 

0.1201 

2.369 

B 

Slope  of  hill 

3 

2.7278 

0.9093 

17.938 

R 

AB  interaction 

4 

0.0659 

0.0165 

0.3255 

AR  interaction 

6 

0.2686 

0.0448 

0.8838 

ABR  interaction 

12 

0.1389 

0.0116 

0.2288 

Error 

179 

9.0475 

0.5056 

From  the  ANOVA,  the  number  of  hills  and  slope  are  significant 
at  the  0.01  level,  but  the  height  of  hills  is  not  significant. 
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b.  Dependent  Speed 


SUMMARY  TABLE  - AN OVA 

1/2 

(X  Transgeneration  Data) 


Source 

DF 

SS 

MS 

F-Ratio 

# of  hill 

2 

8.2224 

4.1112 

79.337 

A 

Height  of  hill 

2 

0.2561 

0.1281 

2.4272 

B 

Slope  of  hill 

3 

2.6384 

0.8795 

16.972 

R 

AB  interaction 

4 

0.06446 

0.01611 

0.311 

AR  interaction 

6 

0.2905 

0.0484 

0.934 

ABR  interaction 

12 

0.1289 

0.0107 

0.2065 

Error 

179 

9.2756 

0.0518 

0.2065 

The  results  of  the  dependent 

speed  case 

are  analogous  to 

those  of  the  constant 

speed 

case . That 

is,  the 

number  of 

hills  and  slope  of  hills  are  significant,  but  the  height  of 
hills  is  not  significant. 

The  mean  and  standard  deviation  of  the  treatment  levels 
are  given  below. 
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SLOPE 


MEAN 


STD  DEV 


CONST 

DEP 

CONST 

DEP 

DOWN 

4.2075 

4.2075 

1.04444 

1.0444 

SIDE 

4.1852 

4.1852 

1.1339 

1.1339 

UP 

3.3148 

3.3333 

1.1947 

1.2131 

BOTH 

4.4259 

4.4259 

1.3818 

1.3818 

# OF  HILL 

MEAN 

CONST 

DEP 

STD 

CONST 

DEV 

DEP 

6 

4.930 

3.038 

1.1871 

0.7166 

12 

4.139 

4.139 

1.039 

1.0388 

18 

3.042 

4.931 

0.7207 

1.1788 

HEIGHT  OF 
HILL 

MEAN 

CONST 

DEP 

STD 

CONST 

DEV 

DEP 

100 

3.875 

3.875 

1.393 

1.393 

200 

4.069 

4.069 

1.191 

1.191 

300 

4.155 

4.169 

1.191 

1.195 

1 
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V.  CONCLUSION  AND  RECOMMENDATION 

Terrain  is  one  of  the  key  factors  in  the  analysis  of 
land  combat.  Thus,  the  representation  of  terrain  in  a land 
combat  model  is  a critical  factor  in  the  model  design.  High 
resolution  combat  models  typically  employ  actual  terrain 
data  to  describe  the  terrain  environment.  This  approach, 
although  realistic  for  the  specific  terrain  being  represented, 
is  very  costly  and  time  consuming  to  implement. 

The  goal  of  the  research  by  Needels  and  this  thesis  is 
to  provide  a parametric  representation  of  terrain  usable 
either  as  input  to  a high  resolution  combat  model  or  as  a 
f "stand-alone"  model  to  investigate  intervisibility  segment 

length  parameters.  In  particular,  the  current  model 
describes  the  terrain  as  a function  of  the  number  of  hill 
masses,  the  height  of  the  hills,  the  position  and  spread 
of  the  hills,  and  the  slope. 

A MBV^I  distribution  function  and  a linear  function 
cutting  the  MBVN  function  are  employed  to  achieve  the 
terrain  representation.  The  terrains  given  in  Figures  12 
through  17  illustrate  the  flexibility  of  the  model.  In 
fact,  with  appropriate  selection  of  the  input  parameters 
(given  in  Figure  9)  any  desired  terrain  may  be  generated. 

The  results  of  the  analysis  described  in  Chapter  IV 
indicate  that  slope  is  highly  significant  in  the  resulting 
values  of  intervisibility  pareuneters . Some  of  the  trends 
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noted  from  the  results  given  in  Chapter  IV  are  described 
below. 

1.  An  up  slope  (cut  from  target  side  to  observer  side) 
resulted  in  the  largest  value  of  percent  intervisibility 
time . 

2.  Increases  in  the  number  emd  height  of  hills  resulted 
in  decreases  in  the  percent  intervisibility  time. 

3.  Speed  dependent  on  terrain  slope  resulted  in  a 
slight  reduction  in  percent  intervisibility  time  from  the 
constant  speed  case. 

The  analysis  '•  f the  hypothetical  data  described  in 
Chapter  IV  serves  to  illustrate  the  potential  utilization 
of  the  model.  The  trends  indicated  from  the  amalysis  are 
appropriate  for  the  selected  terrains  and  conditions,  but 
should  not  be  construed  as  generalized  results. 

An  alternate  method  (tilted  method)  for  generating 
unsymmetric  terrain  outlined  below  is  suggested  as  an  area 
for  future  enhancement  of  the  model.  Given  a MBVN  density 
function  centered  at  Y in  Figure  11,  define  9 as  the  tilted 
angle,  assuming  p = 0 . The  height  of  the  density  function 
is  Z2  at  Y2  eUid  is  given  by 


Z2 


x-y 


MH*EXP{-  I [(  — 


X,  2 


y-y. 


X y 


(38) 


cind 


Z2 

Y2-Y 


tan  [j  - 0] 


(39) 
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FIGURE  IL  TILTED  METHOD 


From  38  eind  39, 


MH*EXP{-  4 [(-^)^  + ( If]}  = (y-y  ) tan(^-0)  (40) 

^ Ojj  Oy 

For  a specified  value  of  x and  y,  the  value  of  ^2 

computed  using  the  golden  section  method.  The  tilted  height, 

ZH,  can  then  be  computed  as  follows; 


ZH 


sin  9 


(41) 


Additional  investigation  of  this  method  is  required  before 
implementation  in  the  model. 

In  conclusion,  a model  has  been  developed  and  exercised 
which  provides  a potentially  valua±ile  tool  for  efficient 
analysis  of  various  terrain  effects  in  the  land  combat 
environment.  The  intervisibility  parameters  output  from 
the  model  can  serve  as  critical  inputs  to  lemd  combat  models 
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for  evaluating  such  things  as  the  contribution  of  mobility, 
agibility,  and  firepower  to  battlefield  survivability. 
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Height;  300, 


FIGURE  13-2.  Terrain:  2,  # of  Hills;  18,  Average  Height; 
300,  Slope;  -Up 
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v^ 


FIGURE  16-2.  Terrain:  3,,#  of  Hill?;  6,  Average  Height; 
300,  Slope;- Side 
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FIGURE  17- 
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